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ABSTRACT 

Subject of this paper are intrinsic ellipticity correlations between galaxies, their statistical 
properties, their observability with future surveys and their interference with weak gravita¬ 
tional lensing measurements. Using an angular momentum-based, quadratic intrinsic align¬ 
ment model we derive correlation functions of the ellipticity components and project them 
to yield the four non-zero angular ellipticity spectra C|(^), Cy,{), Cp(^) and in their 

generalisation to tomographic surveys. For a Euclid-like survey, these spectra would have 
amplitudes smaller than the weak lensing effect on nonlinear structures, but would constitute 
an important systematic. Computing estimation biases for cosmological parameters derived 
from an alignment-contaminated survey suggests biases of -i-5crv„ for the dark energy equation 
of state parameter w, —lOcrn,,, for the matter density Q,„ and -12cro-j for the spectrum normal¬ 
isation (Tg. Intrinsic alignments yield a signal which is easily observable with a survey similar 
to Euclid; While not independent, significances for estimates of each of the four spectra reach 
values of tens of cr if weak lensing and shape noise are considered as noise sources, which 
suggests relative uncertainties on alignment parameters at the percent level. 
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1 INTRODUCTION 

Intrinsic alignments are a fascinating topic: The idea, that proper¬ 
ties and orientation of the luminous stellar component of a galaxy 
depends on the surrounding dark matter distribution is very per¬ 
suasive, but the exact mechanisms and parameters are not well un¬ 
derstood (for recent reviews, please refer to Troxel & Ishak 2014b; 
Joachimi et al. 2015; Kiessling et al. 2015; Kirket al. 2015). On one 
side, tidal shearing is thought to be responsible for the alignment 
of ellipticial galaxies, where the halo potential is distorted by grav¬ 
itational fields of the large-scale stmcture. This affects the stellar 
component as well, leading to a correlation between the brightness 
distribution and the extrenal tidal gravitational field. On the other 
side, and this will be the focus of this paper, tidal torquing aims 
to explain the alignment of spiral galaxies: Here, the stellar disk 
aligns itself with the halo angular momentum direction, which in 
turn has been imprinted on the halo by torquing processes during 
halo formation. This alignment model for spiral galaxies, its ob¬ 
servable signatures and its impact on weak lensing is the topic of 
this paper. 

The mechanism of tidal torquing as a perturbative process in 
structure formation has been quantitatively worked out by Peebles 
( 1969) and was first in detail investigated in numerical simulations 
by White (1984); Warren et al. (1992); Bailin & Steinmetz (2005). 

* e-mail: bjoern.malte.schaefer@uni-heidelberg.de 


Due to the fact that tidal torquing is perturbative (Heavens & Pea¬ 
cock 1988; Catelan & Theuns 1996; Schafer 2009), correlations 
between angular momenta of galaxies can be traced back to corre¬ 
lations in the tidal fields, giving rise to short-ranged correlations, 
with typical correlation lengths of about one Mpc (Crittenden et al. 
2001 ; Catelan & Porciani 2001 ; Schafer & Merkel 2012; Libeskind 
et al. 2013), and correlations between the spin-field and density 
(Lee & Pen 2001). Tidal torquing depends on the environment of 
the haloes, which can be described by extentions to the standard 
torquing picture (Codis et al. 2015). 

Ultimately, tidal torquing leads through the alignments of an¬ 
gular momenta to correlations between the shapes of galaxies, 
which can be quantified by angular correlation functions or angu¬ 
lar spectra, either between the shapes themselves or between the 
galaxy shape and another property of the large-scale structure. 

Evidence for intrinsic alignments and for tidal torquing in 
particular has been found in numerical simulations, from simula¬ 
tions with only dark matter (Altay et al. 2006; Lee et al. 2008; 
Faltenbacher et al. 2009; Aragon-Calvo Sl Yang 2014) to those 
with a baryonic component and star formation (Tenneti et al. 2015, 
2014,?; Schaye et al. 2015; Crain et al. 2015) aiming to provide a 
shape of the luminous commponent. Recently, Dubois et al. (2014); 
Codis et al. (2014) showed that merging processes play an impor¬ 
tant role reorienting the angular momentum direction of haloes, 
which would be in contradiction to the basic idea of predicting an¬ 
gular momentum directions perturbatively. 
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Observations of intrinsic alignments of spiral galaxies are dif¬ 
ficult, and up to now there is not a consistent picture, partially 
because the large number of different quantifiers employed make 
comparisons difficult: While early studies have found correlations 
between density and shapes of galaxies (Lee & Pen 2002; Mandel- 
baum et al. 2006), detection reports have been reported for a range 
of galaxy types in different surveys (Brown et al. 2002; Lee & Er- 
dogdu 2007; Lee & Pen 2007; Paz et al. 2008; Lee 2011; Hirata 
et al. 2007; Siverdet al. 2009; Mandelbaum et al. 2011,2013; Tern- 
pel & Libeskind 2013; Tempel et al. 2013; Joachimi et al. 2013,?; 
Li et al. 2013; Zhang et al. 2013; Singh et al. 2014), whereas (An- 
drae & Jahnke 2011; Heymans et al. 2012, 2013) did not detect an 
alignment signal from these galaxies, but instead one due to ellipti¬ 
cal galaxies. Therefore, one of the motivation of this work was the 
derivation of all observable spectra of the intrinsic ellipticity field, 
and an estimate how well they could potentially be observed with 
the future Euclid mission' (see Laureijs et al. 201 1). 

A second motivation for improving the understanding of in¬ 
trinsic alignments is weak gravitational lensing: Correlations be¬ 
tween galaxy shapes can be due to both effects (Lee & Pen 2000; 
Catelan et al. 2001). If one were to fit a model for a weak cosmic 
shear spectrum to data containing both lensing and intrinsic shape 
correlations, severe parameter estimation biases would be the con¬ 
sequence (Kirk et al. 2010; La.szlo et al. 2012; Kirk et al. 2012; 
Capranico et al. 2013; Valageas 2014). This has been investigated 
in the context of the Dark Energy Survey (DES) and of the weak 
lensing survey of the Euclid-mission, where large biases have been 
estimated to occur due to intrinsic alignments. Shape correlations 
appear on multipoles of ^ 1000 for a survey like Euclid, which 

reaches to redshifts of unity, where the lensing signal is strongest. 

There are a number of ways in which one can deal with in¬ 
trinsic alignment contamination in weak lensing data, from dis¬ 
carding close pairs of galaxies, to nulling (Huterer & White 2005; 
Joachimi & Schneider 2008), and ultimately to provide a param- 
eterised model with enough freedom to enable a fit to the com¬ 
bined spectrum with a minimum of assumptions (Schneider & Bri¬ 
dle 2010). There are already first estimates of the magnitude of this 
effect, for instance from the Mega-Z survey (see Joachimi et al. 
2011 ). 

At the same time are lensing surveys always intrinsic align¬ 
ment surveys and can be used to investigate models of galaxy align¬ 
ment and to determine alignment parameters, if a separation of 
the two effects is possible (Joachimi & Bridle 2010; Kirk et al. 
2010). In fact, gravitational lensing shape correlations and intrin¬ 
sic alignments differ in their statistical and physical properties, for 
instance in S-mode generation (Crittenden et al. 2002), in higher- 
order statistics (Semboloni et al. 2008; Shi et al. 2010; Merkel & 
Schafer 2014), in their appearance in 3d- or tomographic analy¬ 
ses (King & Schneider 2003; Merkel & Schafer 2013) or in their 
cross-correlations to other data sets (Hall & Taylor 2014; Troxel 
& Ishak 2014a). Additionally, there are ways of performing self¬ 
calibrations (Troxel & Ishak 2012b,a; Zhang 2010) and down¬ 
weighting schemes (Heymans & Heavens 2003; King 2005) that 
can control intrinsic alignments as a systematic. Conversely, the 
alignment signal could be enhanced in a statistical way (Joachimi 
& Schneider 2010). 

In summary, we try to answer these questions: Firstly, we aim 
to derive the full set of intrinsic ellipticity spectra resulting from 
an angular-momentum based alignment model (Sect. 3), secondly. 


* http://www.euclid-ec.org/ 


to quantify their observability in a modern tomographic survey 
(Sect. 4), and thirdly, to quantify their interference with parame¬ 
ter estimation from weak lensing (Sect. 5). For the last two points 
we consider the planned Euclid survey, which reaches to a median 
redshift of 0.9. 

The reference cosmological model used is a spatially flat 
wCDM cosmology with Gaussian adiabatic initial perturbations in 
the cold dark matter density and a homogeneous dark energy com¬ 
ponent. The dark energy equation of state-parameter w is constant 
in time. The specific parameter choices are fl„, = 0.25, = 1, 

CTg = 0.8, Dj = 0.04 and Hq = 100 h km/s/Mpc, with h = 0.7, 
and w = -0.9. We adopt the summation convention for implied 
summation over repeated indices. 


2 COSMOLOGY 

2.1 Dark energy cosmologies 

In spatially flat dark energy cosmologies with the matter density 
parameter Q,„, the Hubble function B(a) = d In a/dt is given by 

^+(1-0„,)exp(3 dlna(l -I- w(a))), (1) 

with the dark energy equation of state w(a). We take w(a) to be a 
constant with the value w = -0.9. The relation between comoving 
distance^ and scale factor a is given by 

in units of the Hubble distance ;y// = c///o- 


2.2 CDM structures 

The linear CDM density power spectrum P{k) describes the vari¬ 
ance of the field 6 in the case of homogeneous Gaussian fluctu¬ 
ations, (6(k)d{k'y) = {InySoik - k')P(k), and is given by the 
inflation-motivated ansatz (Bardeen et al. 1986) 

P(k)oc k"‘T^ik), (3) 

with the transfer function T(k). 

The normalisation of the spectrum P(k) is the variance erg on 
the scale R - Mpc//i, 

T r Pdk , 

0-1 = J -^P(k)WfkR) (4) 

where a Fourier transformed spherical top hat filter function, 
W(x) = 3ji(x}/x is used. j{{x) is the spherical Bessel function of 
the first kind of order £ (Arfken & Weber 1995). 

Growth of the density field, 6(x, a) = D+(a)6(x, a = 1), in the 
linear regime |(5| « 1 is described by the growth function D+(a), 
which follows as the solution to the growth equation (Turner & 
White 1997; Linder & Jenkins 2003), 

d-D+{a) 2 - q(a) dD+(a) 3Q.„(a) 

, 9 +-^- = —^D+(a), (5) 

da- a da 2a^ 

where both Cl,„(a) and the deceleration parameter q(a) carry the 
dependence on dark energy. 

The gravitational potential O (normalised by P), which we 
use to predict the two observables gravitational lensing shear and 
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galaxy ellipticity, is related to the growing density field through the 
comoving Poisson equation, 


AO = 


3a„, 5 
2x1 a 


( 6 ) 


such that the gravitational potential grows oc D+(a)la. For gravita¬ 
tional lensing on nonlinear scales we use the model by Smith et al. 
(2003). In contrast, the intrinsic alignment spectra are computed 
from a linear CDM-spectrum, because tidal torquing is a perturba¬ 
tive process and is consistently computed from the initial condi¬ 
tions of structure formation. 


2.3 Weak gravitational lensing 

Weak lensing refers to the correlated change in the shapes of distant 
galaxies by gravitational deflection and distortion by potentials in 
the cosmic large-scale structure (for reviews, see Bartelmann 2010; 
Bartelmann & Schneider 2001), which is one of the prime ways 
to measure cosmological parameters (summarised in Hoekstra & 
Jain 2008; Kilbinger 2014). The two components , and of 
the complex weak lensing shear y, = y+j + iy^j in tomography bin 
i with edges at comoving distances and^,+i provide a mapping 
of the projected second derivatives of the gravitational potential O 
(Takada & Jain 2004; Hu 1999; Takada & White 2004), 


dxW,(x){d^^-dl)o (7) 

X XH 

AxW,{x){2a%)^ (8) 

weighted by the lensing efficiency function 
D+(a) 

w.ix) = G/Ork (9) 

a 

with 

y' _ y 

dx'nix')^^-^. ( 10 ) 

maxO^.r',) T 


Here, nix) - riiz)Azldx = niz)H(z) is the distance distribution of 
the galaxy sample and Xi denote the bin boundaries in terms of 
comoving distance. 

Correlation properties of the two shear components can be de¬ 
scribed by two correlation functions (Kaiser 1992), 

c:,/0) = <7.,,r;y> + <rx,,r'(ii) 

c!,/F) = <y+,,r;j> - <rx,,rxj>- (12) 

Subsequent Fourier-transformation yields the two tomographic 
spectra and C\..iC) (Schneider et al. 2002; Schneider & 

Kilbinger 2007; Fu & Kilbinger 2010), 


Cl,,it) = 

r ede [cl -jie)hm + cl ,^ie)um), 

(13) 

Cl„{t) = 

r ede [c\,-ie)Mm - cl ,^(e)um) ■ 

(14) 

The positive-parity spectrum C^-ff) is the primary lensing 


observable, while the negative-parity spectrum C^g^jit) is zero to 
lowest order in lensing, in contrast to intrinsic alignments, where 
both ellipticity spectra Cgjif) and C’g-it) are nonzero, and in fact 
measurable. 

For everything related to weak lensing, we use the specifica¬ 
tion of a Euclid-like weak lensing survey: median redshift z^ed = 


0.9, density of 40 lensed galaxies per squared arcminute and an el¬ 
lipticity shape noise of cr^ = 0.3. The specific redshift distribution 
niz)Az is parameterised by: 

= .xp(-(i)]fe wich 4 = |r(|)- (15) 

with the parameters yS = 3/2 and za = Zmcd/ V2, and tomographic 
redshift bins are chosen to contain an equal fraction of the total 
number of galaxies. Likewise, the bin boundaries^, for weak lens¬ 
ing tomography, and later for the ellipticity correlations, are cho¬ 
sen to contain equal fractions of the total galaxies, which keeps the 
shape noise in each bin constant. 


3 INTRINSIC ALIGNMENTS 
3.1 Angular momenta of galaxies 

The physical picture for angular-momentum induced alignments in 
spiral galaxies is the following: The host CDM halo picks up angu¬ 
lar momentum through tidal torquing (Peebles 1969; White 1984) 
prior to gravitational collapse. In this process, both baryonic and 
dark matter have equal amounts of specific angular momentum, 
and when at a much later stage a galactic disk is formed, its sym¬ 
metry axis should be aligned with the host halo angular momen¬ 
tum axis, if neither dissipative processes nor merging activity have 
have changed the angular momentum or the disk orientation. Com¬ 
monly, one uses the alignment model by Lee & Pen (2000, 2001) 
for setting up a distribution p(L|Oayj)dL for the angular momentum 
direction L conditional on the tidal shear field 5^^® through the 
covariance matrix, 

{L„Lp) = - I —^— a (16) 

where is the unit-normalised traceless tidal shear. The param¬ 
eter a interpolates between isotropic and tightly coupled angular 
momentum directions and is measured to be a 1/4 in numeri¬ 
cal simulations for galaxy-sized haloes. It effectively describes the 
amount of misalignment between the gravitational shear tensor and 
the moment of inertia of the halo prior to collapse (Porciani et al. 
2002a,b). We emphasise that our model links angular momentum 
and ultimately galaxy shape to the squared tidal shear in an av¬ 
eraged way without distinguishing between environments whose 
dynamical properties can have an influence on the spin statistics 
(see, for instance, Kashikawa & Okamura 1992; Cuesta et al. 2008; 
Aragon-Calvo et al. 2007; Codis et al. 2015). Assuming perturba¬ 
tion theory allows to link the angular momentum direction to the 
initial conditions of structure formation: This assumption has been 
tested by Porciani et al. (2002a,b); Lee & Pen (2008) in simula¬ 
tions, who found that the angular momentum direction can be pre¬ 
dicted well, in contrast to the angular momentum magnitude. 


3.2 Angular momentum alignment 

The idea of angular-momentum induced alignments (so-called 
quadratic alignments) is an alignment of the symmetry axis of 
the galactic disk with the host halo angular momentum direction 
(Heavens et al. 2000; Crittenden et al. 2001, 2002; Mackey et al. 
2002). The ellipticity e = e+ + icx measured by an observer then 
depends on the direction L = L/L of the angular momentum vector 
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L relative to the line of sight: 


( 22 ) 


Ll-Ll 


1 +L? 


and 


Ex 


= 2a 


LjxLjy 

i + U’ 


(17) 


if the line of sight is taken to be the z-direction. The modulus \e\ of 
the ellipticity depends on L according to 


\e\ = a 


l-L? 

l+Lf 


( 18 ) 


by using the normalisation = 1 of the angular momentum direc¬ 
tion. a is a phenomenological parameter which weakens the depen¬ 
dence of e on L if taken to be smaller than one, and in this work we 
adopt the conservative value of cr = 1/2. Pictorially, the parameter 
describes the effect of a galactic disk of finite thickness, where the 
change in ellipticity with inclination angle is weaker compared to 
an idealised infinitely thin disk. Crittenden et al. (2001) have com¬ 
puted the mean value of the complex ellipticity e for a given tidal 
shear averaged over angular momentum magnitudes and directions, 
and arrive at 

e(®) " Y 2 j “ ®r«®<»r “ 2i0.,„0„,). (19) 


Similarly, the absolute value e, = |e| = + Cx of th® complex 

ellipticity e = i€x can be related to the tidal shear through 

3 /X 

es = ^aa2j®ztt®az, (20) 

a 

using the normalisation condition = 1. It should be em¬ 

phasised that the scalar ellipticity e, is not statistically independent 
from the complex ellipticity e - e++ ifx and that it will in fact con¬ 
tain less information than the complex ellipticity due to removal of 
the ellipse’s phase angle, i.e. its orientation, while only keeping in¬ 
formation about the axis ratio. This type of information is as well 
a source of lensing information as shown by Heavens et al. (2013); 
Alsing et al. (2014) 


3.3 Ellipticity correlations 

Neighbouring galaxies have formed from similar initial conditions 
and show ellipticity correlations through correlations in the tidal 
shear they have been exposed to while building up their angular 
momenta. It is important to emphasise that alignment processes are 
thought to be local and do not arise through interaction between 
the galaxies: The are generated only through correlations in the 
aligning tidal field 3^^®. Therefore, correlations in e must be ex¬ 
pressed in terms of correlations of squares of the tidal shear ®. 
This proceeds in two steps. Firstly, the tidal shear field correlations 
Cti^yd(r) = (®Q.^(a:)®ya(x')), r = \x - x'\, can be computed to be 

Cafiysir) {,^ap^y6 ^ay^pS ^a6^p,y )(( 2 ('') + 

darpSys + 5 perm.)^ 3 (r)-l- (21) 

r^rpryKs^^ir) 

from which the correlations Capys(r) of the traceless shear follow by 
subtraction of the trace ®Q.a/3d„^ from the tidal shear <!>ap. The unit 
vector is chosen to have the entries (sin a, 0 , cos a) by suitable 
rotation of the coordinate frame. Finally, correlations of the squared 
tidal shear components are related to Capysir) by virtue of 

(®a(x)®b(x) ®c(x')®;,(x')) = [CacCbd + CadCbc) , 


where the capital letters denote the index pairs of the tidal shear 
tensors ®„^(x). 

We generalise the derivation presented by Crittenden et al. 
( 2001 ) to include all possible correlations between the ellipticity 
components e+, 6 x and e^. They follow from contraction of the ex¬ 
pressions (19) and ( 20 ) while substituting r„ for our orientation 
of the coordinate frame and while making use of the Kronecker- 
6 symbol. 


(e+<> 

1 ( aa f 

= Y 4 \ 14^ I cos(4(i') - 1 - B++ cos( 2 a) -l- C++ 

) (23) 

(fxO 

1 / aa 

(24) 

(e+e,) 

1 1 aa f 

= Y 24 1 14 ^ ) cos(2a) - 1 - C+,) 

(25) 

(gg> 

1 ( aa f 

= 108 1 14^ 1 ‘^°s( 2 Q') + Cj.,) 

(26) 

with the abbreviations 



(-4 + 6(k3 + + llil 

(27) 

Byy = 

- 32^3('4 - 28/'^ -H 

(28) 

Cyy = 

+ (26('3 -1- 58('2)('4 -H 155^^ -n 472^:2^3 + 336^2^ 

(29) 

for the correlation function (e+e(). 


zlxx = 

0 

(30) 

Bxx = 

(C3 - 3/-2)C4 - 5(1 - 

(31) 

Cxx — 

(-(-3 - 5(2)(a - ISCf - 59(2(i - A2(l 

(32) 

for the function {e^e'y.) with Axx being zero. These results have been 
derived by Crittenden et al. (2001), and are listed for completeness. 
In addition, we supply the new results 

-^-FS “ 

-5(1 + 18(('3 -h ^2)C4 + 51(1 

(33) 

— 

4/-2 -H (32(2 + AS(2)(4 + m(l + A0S(2(3 

(34) 

C+s = 

1(1 + (-50^3 - 66f2)('4 - 233(1 - 408(-2('3 

(35) 

for the correlation (e+ej) and lastly 


^ss — 

9(1 + (5A(2 + 5A(2)(4 + 153(1 

(36) 

Bss = 

60(1 + (480^3 + 288 ( 2 X 4 + 1284^32 -H 1800('2C3 

(37) 


C„ = + (490^3 + 426^2)^4 + ^901 + 3864/'2C3 2448f| 

for {Psd)- We checked that the remaining correlation functions 
(e+ex) and (exfi) are in fact zero. The function f„(r), 

(- 1 V r dk 

J r-dAkr)P{k) (38) 

denotes weighted moments of the CDM-spectrum P{k). Here, we 
impose a Gaussian smoothing on the CDM-spectrum P(k), 

P(A:)^P(fc)exp|-^J, (39) 

where the smoothing scale R is linked to the mass scale of the ob¬ 
jects under consideration through M = AnlhQ-^pcruR^. Specifically, 
we set the mass scale to M = IO'^MqIIi, which is typical for a spi¬ 
ral galaxy. Because spin generation in a halo is taking place prior 
to gravitational collapse, the smoothing scale uses the cosmic value 
DmPcrit for the density. We work with a Gaussian smoothing for nu¬ 
merical purposes. 

In the very last step. Limber projection of the 3d-correlation 
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Figure 1. Ellipticity correlation functions C' ^(8) (solid lines) and C‘ -{d) 
(dashed lines) for 4-bin-tomography with a Euclid-like survey. 


functions (ee')(r,a) with the galaxy distribution nix) in the corre¬ 
sponding tomography bin i yields the angular ellipticity correlation 
functions {ee')(d) for all four non-zero combinations of ellipticity 
modes. The angular correlation functions depend on a single bin 
index, i.e. there are no cross correlation of the ellipticity field be¬ 
tween different bins, because the correlation length of the ellipticity 
field is much shorter than the bin width: Typical numbers would be 
a few Mpc for the ellipticity field in comparison to hundreds of 
Mpc for Euclid’s tomography bin size. 


3.4 Ellipticity correlation function 

The two components of ellipticity can be combined into a complex 
ellipticity e = -t-ifx, such that in analogy to weak lensing or to the 
polarisation of the cosmic microwave background two correlation 
functions C' .(6) and Cf -^6) can be defined: 

Cip) = <e+,,e;,)-f <ex,,ei,,), (40) 

Cj(e) = <e+,,e;,> - (41) 

These correlation functions can be complemented by those involv¬ 
ing the modulus e, of the ellipticity, a cross correlation C^,(6) and 
an auto-correlation C| .(0): 

Cim = <£+,,<,■>, (42) 

Cim = <£„■<,-). (43) 

In contrast to lensing, in which correlations between different to¬ 
mography bins are natural because the light from galaxies inside 
a distant bin needs to transverse foreground bins which gives rise 
to a correlated lensing signal, the ellipticity correlations only occur 
locally, due to the short-rangedness of tidal fields. 

Fig. 1 shows the two correlation functions C‘1 .(9) and Cf .(9) 
for a 4-bin tomographic survey. Clearly, intrinsic ellipticity corre¬ 
lations are short-ranged and only present on scales below a few 
arcminutes. In both cases the ellipticity field shows correlations of 
the order 10“^ on small scales, which is smaller than that caused by 
gravitational lensing. Only on very small scales the two effects will 
become comparable. 

The correlation functions C^.(9) and .(6) are depicted in 
Fig. 2. Again, sizeable amplitudes are only present below the scale 


Figure 2. Ellipticity correlation functions C^,(6) (solid lines) and CJ^(6) 
(dashed line) for 4-bin tomography with a Euclid-like survey. The scaling 
has been chosen to make this figure directly comparable to Fig. 1 . 

of a few arcminutes, and fall short in comparison to those of gravi¬ 
tational lensing. 

3.5 Angular ellipticity spectra 


In complete analogy to weak lensing, in total four nonzero elliptic¬ 
ity spectra can be derived from the ellipticity correlation functions. 
The first pair is 


cyt) = kJ 

9d9 [cy9)Jom + cyey^m), 

(44) 

cyt) = ^ j 

9A9 [cy9)Jom - cy9)j^m) , 

(45) 

from the correlation functions C^.(0) and Ci-iO) and the second 
pair 

Cyt) = 2;r 

I' M9cy9)j2(m, 

(46) 

Cy£) = 2;r 

1' 9d9cy9)Jom, 

(47) 


while all other combinations vanish due to statistical parity invari¬ 
ance. 

The spectra -(£) and Cg are shown in Fig. 3 alongside 
a grey band which shows the amplitude of the weak lensing spec¬ 
trum Cg from linear and from nonlinear structures, computed 
for a 4-bin survey. At low multipoles both spectra are constant and 
equally large, only at multipoles above a few hundred the spectrum 
Cgff) dominates over the spectrum Cg({) by about an order of mag¬ 
nitude. The shape of the spectrum depends on the redshift interval 
of the tomography bin because correlations of the same physical 
correlation length appear through the choice of redshift bin on dif¬ 
ferent angular scales. Ultimately, the spectra drop of rapidly due 
to the exponential smoothing of the CDM-spectrum. It is interest¬ 
ing to not that, with our conservative choice for the misalignment 
parameter a and the disk thickness a intrinsic alignments provide 
a significant contribution to the £-mode weak lensing spectrum, 
which will be the focus of Sect. 5. 

Finally, the spectra ■{() and C| ■{£) are given in Fig. 4, again 
in comparison to the weak lensing spectrum C'^.j(£) for a 4-bin 
survey. They tend to have smaller amplitudes compared to 
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Figure 3. Angular ellipticity spectra .{€) (solid lines) and .{£) (dashed 
lines) in comparison to weak lensing spectra CZ ■ (0 (shaded area) for linear 
(light grey) and nonlinear CDM-spectra (dark grey), for 4 tomography bins. 
The colour gradient ranges from blue for the low-redshift bins to green for 
the high-redshift bins. 



10 30 100 300 1000 3000 

multipole £ 


Figure 4. Angular ellipticity spectra C^j(0 (solid lines) and C| -(0 (dashed 
lines) in comparison to weak lensing spectra CZ ■ (0 (shaded area) for linear 

£.,lj 

(light grey) and nonlinear structure formation (dark grey) and for 4 tomog¬ 
raphy bins. The colour gradient changes from blue for the low-redshift bins 
to red for the high-redshift bins. 

or and in particular has much smaller amplitudes at 

low multipoles. 


4 OBSERVABILITY 

A natural question concerns the significance at which intrinsic 
alignments can be observed in future, possibly tomographic sur¬ 
veys. For this purpose, we compute a forecast for the signal to 
noise-ratio which can be reached with a Euclid-like survey un¬ 
der the assumption of Gaussian statistics and a perfect separation 
between weak lensing and intrinsic alignment correlation. Noise 
sources considered are shape noise in the ellipticity measurement 
and cosmic variance. In a previous paper (Capranico et al. 2013) 



Figure 5. Signal amplitude for observing the ellipticity spectra ■{{) 
(solid lines) and CJ -(t) (dashed lines) as a function of multipole { and the 
total number of redshift bins into which the survey is divided. The survey 
was assumed to cover the full sky. The signal amplitude increases with the 
number of tomography bins. 

we have shown at least for a non-tomographic survey, that using 
priors on cosmological parameters from e.g. observations of the 
cosmic microwave background and from baryon acoustic oscilla¬ 
tions allow a sufficiently accurate prediction of the weak lensing 
contribution of the spectrum such that it can be subtracted. 

We estimate the signal to noise-ratio using the Fisher- 
formalism for determining the error on the amplitude of the signal 
and compute the signal to noise-ratio as the inverse relative error 
on the unknown normalisation of the signal. Specifically, the signal 
strength E follows then from 

Y? = ^yp.l+\)\x(c-'S C-'S) (48) 

t 

with the signal covariance S,/ and the total covariance C,;, which 
contains in addition to the intrinsic ellipticity spectra the contribu¬ 
tion from weak gravitational lensing and ellipticity shape noise. In 
this way the cosmic variance, which is mostly generated by gravi¬ 
tational lensing, is properly taken care of. The sum is extended over 
all multipoles up to f = 3000 to capture most of the signal before 
the shape noise contribtion dominates. 

Specifically, we consider the case of separately measuring the 
positive parity spectra .{() and C| .(f) first. As shown by Fig. 5, 
the signal increases with increasing multipole E until it levels off 
at ^ 1000 when the shape noise of the ellipticity becomes domi¬ 

nant. Subdivision of the galaxy samples in tomographic bins is able 
to boost the signal significantly from Ys = 6cr (3.4 it) to Ys = 60cr 
(34cr) in the case of CJff) and from Ye = 25cr (14cr) to close to 
Ye = 200(T (114cr) in the case of C‘^(E). The number in parantheses 
correspond to Euclid’s sky coverage of /sky = 0.33, which reduce 
the full sky significances by about 30%, and we considered tomog¬ 
raphy with up to 9 bins in redshift. 

Similar results can be obtained for a combined measurement 
of all three positive-parity spectra C‘^.{E), ^{E) and C^/f)- Due 

to the fact that e, is derived from the two modes 6+ and these 
measurements are not independent and their covariance needs to be 
incorporated into the estimate. Eig. 6 shows that this measurement 
is able to yield Ysec - 200cr (140cr) when using tomography. The 
negative-parity spectrum ■{{), which is predicted to be zero for 
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Figure 6. Signal amplitude for observing the ellipticity spectra Cg ^{t) (solid 
lines) and all three positive-parity spectra C| -{i) and com¬ 

bined (dashed lines) as a function of multipole i and number of tomographic 
bins. The sky coverage is set to one. The signal increases with increasing 
number of tomography bins. 


gravitational lensing, can be easily measured with a significance of 
about Sg = lOOcr (70 it). Clearly, the lack of cosmic variance due to 
gravitational lensing makes up for the intrinsically smaller signal. 

In summary, we find that ellipticity alignments due to corre¬ 
lated angular momenta should yield a signal which is easily mea- 
sureable in future tomographic surveys. As the signal is propor¬ 
tional to the alignment parameter a and the disk thickness param¬ 
eter a, one can expect errors on the percent-level on these two pa¬ 
rameters of the quadratic alignment model. It should be kept in 
mind, however, that the product of the parameters determines the 
amplitude of the intrinsic ellipticity spechum and that they can not 
be measured separately. 

From the cosmological parameter set, only fl,„, h, n, and to 
a lesser extend w influence the shape of the ellipticity correlation 
function, while erg is not relevant in the quadratic model: This can 
be seen in eqns. (19) and (20), where the ellipticity is linked to the 
traceless tidal shear which does not depend on the amplitude of 
fluctuations in the tidal shear field as a consequence of the angular 
momentum model. Secondly, all terms in the ellipticity correlation 
functions depend on ratios between -functions such that the mag¬ 
nitude of is canceled. Q.,,,, h and n, determine the shape of the 
CDM-spectrum and influence therefore the shape of the ellipticity 
correlations through the ^„-functions. Additionally, the shape of the 
angular correlations is affected by the conversion between redshift 
and comoving distance in the Limher-projection, which depends on 
n„, and w. 

We conclude that intrinsic alignments are by far the most im¬ 
portant secondary effect in weak lensing surveys even at high sur¬ 
vey depth, and surpass in particular in 6 -mode generation of other 
effects such as Bom-corrections or clustering. Baryonic effects on 
the matter spectrum are of similiar order compared to intrinsic 
alignments. 


5 PARAMETER ESTIMATION 

Weak lensing shear operates on the complex ellipticity e = e+ + iex 
of a galaxy through the mapping e —> e + y in the limit of small 


shears. Therefore, correlations of the observed ellipticities will con¬ 
tain the intrinsic ellipticity correlation (ee'), the correlation of the 
lensing shear ( 77 ') and possibly cross-correlations between intrin¬ 
sic shapes and gravitational shear, (ye' +y'e), although the last term 
is zero in the case of the quadratic alignment model and Gaussian 
fluctuations. Therefore, it is possible to add lensing and intrinsic 
alignments as independent contributions to the total ellipticity spec¬ 
trum. 


5.1 Statistical errors 

The standard method of computing forecasts of statistical errors 
is the Fisher-matrix technique (Tegmark et al. 1997; Bassett et al. 
2009,2011 ), in which the cross-correlation between different tomo¬ 
graphic bins can be correctly incorporated. We derive forecasts on 
cosmological parameters from the positive-parity shear spectmm 
Cg -jXt) and use the expression 

^ y (2^ + 1) tr (d^ In C a. In c) (49) 

f 

for the Fisher-matrix, where the scaling oc lowers the signal 
amplitude due to incomplete sky coverage of the survey. In the case 
of tomographic bins which contain equal fractions of the galaxy 
sample, the noise to be added to the covariance is diagonal, 

C, 7 (^) = + ^ 0 -%. (50) 

From the Fisher-matrix it is straightforward to define condi¬ 
tional errors ^^d marginal errors = (F“*)^^ on 

cosmological parameters, where we use the set Q^, erg, h, and 
w, while imposing spatial flatness. Commonly, we derive statistical 
errors hy summing over all multipoles until £ = 3000, and vary the 
number of tomography bins between «bin = 2 and 5, which results 
in a lensing signal of close to Z = lOOOcr and marginalised statisti¬ 
cal errors on the cosmological parameters on the percent level. 


5.2 Systematic errors 

The presence of intrinsic alignments in weak lensing data would 
lead to biases in the estimation of cosmological parameters if they 
are uncorrected. For this case we consider biases that would arise if 
an alignment contribution , (0 is added to the weak lensing spec¬ 
tra Cg ■■{£) for equal bin indices, while i ■=£ j, is unchanged. 

A fit of Cg .j(£) to the combined .j(£) + (no summation 

over i implied) would then give rise to a parameter estimation bias 
d/j. These estimation biases can be computed from the hue model 
C,_,y = C^-.j + C‘^.6ij and the false, incomplete model Cfjj = 
by solving the linear equation 

2 = S^ = ^(G^')^,,a„ (51) 

V V 

with the vector a^, 

fl,. = y ^^^tr[ajc^lnC/(id-C^‘C,)], (52) 

t ^ 

and the matrix G^y, 

G,y = Yj^-^tv[c-/dlyCf[c-f^C,-iA)] 

i ^ 

- (53) 
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Figure 7. Marginalised likelihood contours (blue-green) and conditional 
likelihood contours (red-blue) at Icr-confidence for the parameter pair 
and (Tg, for Wbin = 2 (smallest) up to 6 (largest), in comparison to parameter 
estimation biases due to intrinsic alignments with a = 0.25, all for fmax = 
3000. 


Figure 8. Marginalised (blue-green) and conditional (red-blue) Icr-contours 
in the n„,-w-plane, along with parameter estimation biases. The colour in¬ 
dicates the number of tomography bins used in the analysis, ranging from 
tibin = 2 (smallest) to 6 (largest), for fmax = 3000. 


where id refers to the unit matrix in tibin dimensions. The formal¬ 
ism employed here (Schafer & Heisenberg 201 2) is a generalisation 
of the bias estimation formalism proposed by Cabre et al. (2007), 
Taburet et al. (2009), Amara & Refregier (2008) and March et al. 
(20 11 ) to tomographic data. It has been shown by comparison to the 
results from Monte-Carlo Markov-chains that it predicts estimation 
biases well, at least in the case of small systematics (Taburet et al. 
2010 ). 

Fig. 7 shows parameter estimation biases in comparison to Icr- 
likelihood contours on Q,„ and erg, both marginalised over the entire 
parameter set and assuming that all other parameters are perfectly 
known. Intrinsic alignments bias both erg and Q„, high by a few 
percent, which should be see in comparison to the very good sta¬ 
tistical precision. The physical reason for biasing is the addition 
of power on the spectrum by intrinsic alignments, which needs to 
be compensated by a fit through increasing the amplitude of the 
signal, leading to larger values of both parameters. Another inter¬ 
esting trend is a growing bias in the cosmological parameters when 
the number of tomography bins is increased. 

Fig. 8 repeats the above analysis for the cosmological param¬ 
eters Cl,„ and w by showing the biases in relation to marginalised 
and conditional likelihood contours. In particular w is measured 
more negative than the fiducial model, hence one could potentially 
mistake a dark energy model for ACDM, again because the intrin¬ 
sic alignments increase the amplitude of the ellipticity spectrum. 
Again, one finds increasing biases if the lensing survey is divided 
into more redshift bins. 

These results are presented by Figs. 9 and 10 in a more quan¬ 
titative way, by illustrating the dependence of the ratio for 

the conditional and marginalised statistical error, respectively. 

Naturally, the biases are much larger than dp/cr^,^ be¬ 

cause statistical uncertainties are smaller in models with fewer pa¬ 
rameters. But in general, the entire ACDM- or wCDM-parameter 
set is affected by intrinsic alignments. Which parameters suffer 
most depends on the number of tomography bins being chosen. 
With 5 bins and = 3000 it is Q,„, erg and that are most 
strongly affected, followed by the dark energy equation of state w. 

Comparing the bias to either the marginalised or conditional 



number of tomography bins 


Figure 9. Systematic error <5^ in units of the conditional statistical error cr^ ^ 
in a wCDM parameter set consisting of Cl,„, erg, h, n, and w, as a function 
of tomography bins nbjn and for fn,ax = 3000. The shaded area indicates the 
3cr-interval. 

error in projections, as done by Figs. 7 and 8, does not give a fair 
representation of biasing relative to the magnitude of the statistical 
errors and their degeneracies. For that reason, we define the figure 
of bias ^ as a quadratic form, 

(54) 

liv 

which compares the systematical error d^ with the statistical one 
and reflects the orientation of the vector d^ with the shape and size 
of the Fisher-ellipse in the full parameter space. In the limit of ab¬ 
sent degeneracies, q corresponds to the quadratic sum of the bias 
components normalised by the statistical errors on individual pa¬ 
rameters. At the same time, exp(-(7“/2) is the evalutated likelihood 
at the position d^ relative to the fiducial model. 

We think that in this way ^ is a convenient way of quantifying 
the total amount of estimation bias relative to the total statistical 
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number of tomography bins 


Figure 10. Systematic error 5^ in units of the marginalised statistical error 
in a h’CDM parameter set consisting of n„,, erg, h, and w, as a 
function tomography bins nt|j„ and = 3000. The shaded area indicates 
the Icr-interval. 

uncertainty of a measurement, but has the disadvantage of not con¬ 
taining any information about the biasing direction as eqn. (54) is a 
positive definite quadratic form. Values of q are typically == 20 for 
{ = 3000 with a weak variation with the number of bins, so rather 
in accordance with the ratio d^/cr^.c, where the bias was expressed 
in units of the conditional statistical error, than with with 

the marginalised statistical error. 


6 SUMMARY 

Subject of this paper are angular spectra of angular-momentum in¬ 
duced, quadratic alignments of galaxies as they would appear in 
tomographic weak lensing surveys, in particular in that of a Euclid- 
like survey. Starting from tidal torquing as a model for generat¬ 
ing galaxy spin we derive the statistics of galaxy ellipticities for a 
given tidal shear and compute the ellipticity correlation function re¬ 
sulting from the tidal shear correlations. Projection and subsequent 
Fourier-transform yield the angular ellipticity correlation function 
and ultimately the angular ellipticity spectrum, whose statistical 
properties we analyse in detail. 

Very important simplifications in our analysis include the pre¬ 
diction of the angular momentum direction by linear tidal torquing 
and the description of the angular momentum field to be derived 
from the tidal shear field in a Gaussian, local random process. Fur¬ 
thermore, we assume that the baryonic component aligns itself per¬ 
fectly with the angular momentum axis of the host halo: This, in 
particular, is a serious simplification which we need to impose in 
our analytical work, such that the derived signal to noise-ratios 
should be considered as upper limits. 

(i) Ellipticity correlation functions show that the ellipticity field 
is only correlated on scales smaller than a few areminutes due to 
angular-momentum induced alignments, and the ellipticity corre¬ 
lation quickly vanishes on scales larger than ten areminuts for all 
components of the ellipticity. 

(ii) Angular ellipticity spectra follow from ellipticity correla¬ 

tion functions through Fourier-transform. We isolate the four non¬ 
zero spectra: the £-mode spectrum the associated S-mode 

spectrum the spectrum CJ(^) of the ellipticity modulus and 


finally the cross-spectrum C^(£) of the £-mode and the elliptic¬ 
ity modulus. These four spectra correspond formally to the four 
non-zero temperature and polarisation spectra of the cosmic mi¬ 
crowave background. Ellipticity spectra are dominated by weak 
gravitational lensing on small multipoles, but have similar ampli¬ 
tudes on multipoles of a few hundred. 

(iii) The significance at which the spectra can be measured was 
estimated assuming Gaussianity of the ellipticity and the weak lens¬ 
ing fields. Euclid turns out to be an excellent mission for investigat¬ 
ing intrinsic alignments: Due to the fact that intrinsic alignments are 
effectively uncorrelated over large distances, each redshift bin pro¬ 
vides a statistically independent estimate of the correlation func¬ 
tion, only with larger Poissonian shape noise. Individual angular 
ellipticity spectra can be measured at significances for the spec¬ 
trum CgCf) of ~ lOcr for a non-tomographic survey increasing up 
to values in excess of lOOtr if the survey is divided into 9 redshift 
bins. Similarly, the ellipticity modulus correlation can be measured 
with ~ 6(T, and this number can be boosted to ~ fiOcr by doing 9- 
bin-tomography. Corresponding numbers for the 6-mode spectrum 
Cp(€) are ~ 20cr and ~ 90cr, respectively. To put these significances 
in context, we would like to emphasise that tomographic weak lens¬ 
ing generates a signal close to lOOOcr. 

(iv) While the four angular ellipticity spectra 

C^it) and Cg({) are not statistically independent because they are 
derived from only two components of the ellipticity field, the com¬ 
bination of all positive-parity spectra would yield a total signal of 
almost 200 it using 9 tomographic bins and extending the estimates 
to multipoles of f = 3000. Due to the fact that signal to noise-ratios 
are inversely proportional to the relative error on the normalisation 
of the signal, we could expect percent-level errors on the alignment 
parameters a and a. 

(v) Estimates of cosmological parameters from weak gravita¬ 
tional lensing are severely biased by intrinsic alignments. While 
the absolute values of the biases do not change much when increas¬ 
ing the number of tomography bins, the bias in units of the statis¬ 
tical error increases significantly: We find biases of -t-Scr^. for the 
dark energy equation of state parameter w, -20<Tfi„, for the matter 
density Q,„ and -lIcTa-g for the spectrum normalisation erg. These 
biases have important implications for measurements of parameters 
as well as for the selection of cosmological models, as one might 
prefer ACDM over a true dark energy model. 

In future work we aim at supplementing the angular- 
momentum based alignment model, which is thought to be applica¬ 
ble to spiral galaxies, with a tidal shearing model for treating ellipti¬ 
cal galaxies (Hirata & Seljak 2004, 2010; Blazek et al. 2012, 2011, 
2015), and combine both models for a realistic morphological mix 
of galaxy types. This would have important consequences: Due to 
the fact that there should be no intrinsic cross-correlation between 
the shapes of spiral and elliptical galaxies for Gaussian random 
fields one can expect smaller levels of intrinsic alignments. If q is 
the fraction of spiral galaxies in a sample, the total ellipticity corre¬ 
lation function would be (ee') = q^{ese[) -l- (1 - q)~{ei.e'P with a cor¬ 
relation function (fjej) of the shapes of spiral galaxies and (fee') of 
elliptical galaxies, respectively. Both factors q^ and (1 - ^)^ would 
be smaller than one, thus reducing the amount of ellipticity corre¬ 
lation relative to the assumption of a single alignment model. Con¬ 
versely, intrinsic alignments have been observationally confirmed 
for elliptical galaxies (for instance by CFHTLenS, Kitching et al. 
2014; Fu et al. 2014), and those have been demonstrated to have an 
impact on parameter estimation, in particular in models with mod¬ 
ified gravity (Laszlo et al. 2012; Dos.sett et al. 2015). The impact 
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of alignments of elliptical galaxies on parameter masurement from 
weak lensing has been estimated by Joachimi et al. (2011) for the 
Mega-Z sample, and biases were found to be significant, even for a 
rather small data set. 

A second interesting topic are cross-correlations between the 
shapes of elliptical galaxies and weak gravitational lensing, which 
are absent in the case of spiral galaxies, due to the same argument 
as before: Lensing, as the shapes of elliptical galaxies, are linear in 
the tidal shear. These cross-correlations are negative and shape the 
ellipticity spectrum in a more complicated way with consequences 
for biases on cosmological parameters. Furthermore, both types 
of intrinsic alignments are local and should not show correlations 
across different tomography bins, and only quadratic alignments 
should be able to excite S-modes at lowest order. All these proper¬ 
ties would be important validations of the basic ideas behind galaxy 
alignments. 
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